Gluon-decay cascade in AA collisions at the LHC and saturation 
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We extract the energy-dependence of the gluon-jet decay cascade from e e~ anni- 
hilation data and show that this energy-dependence has power-low behavior for the 
average transverse momentum of jet bigger than 1 GeV. We show that the observed 
different power-law energy-dependence of hadron multiplicity in AA compared to pp 
collisions at the LIfC is due to the enhanced gluon-decay effect before hadronization. 
This effect is more important for AA collisions where the saturation scale is larger 
than 1 GeV. Here we confront our published predictions in arXiv: 1102.2385 for the 
rapidity distribution of the charged hadron production with the recently released 
data from the CMS collaboration in AA collisions and show that our predictions is 
in perfect agreement with the data. 



1 Introduction 

One of the most unexpected new feature of the LHC data has been the very different 
power-law energy behavior of the charged hadron multiplicities in AA compared to pp 
collisions [TJ [2] ■ There have been several approaches to address the above problem [31 
m [51 [B] . This is closely related to the open problem of entropy production at the early 
stage of heavy-ion collisions (Mini Bangs) . 

In the following, we provide a simple explanation of this observation and confront 
our published predictions [3] with the recently released data from the LHC [2j for the 
rapidity distribution of the charged hadron production. 



2 



A. H. Rezaeian Low-x Meeting 2011 



2 Main formalism and predictions 

In the Color Glass Condensate (CGC)/saturation based approaches (for a recent review 
see Ref. [7]), the hadron production may be divided in two stages: production of giuons 
and subsequently the decay of gluon-jet (or mini-jet) into hadrons. Therefore, the multi- 
plicity of the produced hadrons at pseudo-rapidity rj can be calculated as a convolution 
of these two stages (see Fig. [1]): 
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where '^j^ -a — gives the gluon jet production yield at rapidity y (in pp or AA collisions) 



dyd'-'pT 

computable in the kr factorization scheme [3J [HI [TO] • The second term NGiuon ^^^^ 
average multiplicity of hadrons in the gluon jet with a jet energy Ej^t computable in 
the Modified Leading Logarithmic Approximation (MLL A) scheme [8j [11] . The symbol 
® indicates a convolution, that is, integrals over variables with possible weight factors 
included [3]. 

The greatest simplicity of the CGC approach is the fact that all complexity of infinite- 
body problem at very high energy (or small x) shall be reduced into a one scale problem, 
with a hard saturation scale Qs as the only dimensional relevant scale at which nonlinear 
giuons recombination effects start to become important. In this framework, the secondary 
hadrons are originated from the decay of gluon mini-jets with the transverse momentum 
equal to the saturation scale Qs [U [S] [9]. Then Eq. ([2]) up to a possible logarithmic 
correction, can be readily obtained by a dimensionality argument where Cs is the effective 
area of interaction. 

Notice that the /ct factorization has been proven at the leading log(l/a;) approxima- 
tion for scatterings of a dilute system on a dense one (such as proton-nucleus collisions) 
and includes BFKL type gluon emissions (with gluon fusion effects) between the projec- 
tile and target, and also gluon radiations from the produced giuons [TO]. However, the 
other contribution of the gluon decay, before hadronization, stems from the kinematic 
region outside the BFKL emission regime where both emitted giuons are collinear to 
the emitter [3]. In this kinematic region, the angle between the gluon (quark) and the 
decay gluon is small and the main contribution of gluon-decay has an opposite angular 
ordering to the BFKL type gluon emissions. The later kinematic region has been tra- 
ditionally studied within a different resummation scheme, the so-called MLLA where it 
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Figure 1: The corresponding diagram to the factorization given in Eq. (H]). In the kr 
factorization, the gluon production can be calculated as a convolution of two unintegrated 
gluon density (j)G of the projectile and target ^lOj . 

contains systematically next-to-leading logarithmic corrections and incorporates single 
and double- logarithmic effects in the development of parton cascades [11] . Notice that 
the gluon decay effect incorporated within these two resummation schemes have different 
kinematics [51 [TU1[TT] . 

It is well-known that the gluon decay probability can be factorized from the rest 
of cross-section in e+e^ — > qqg reaction 11 . This is in a sense, the essence of the 
factorization given in Eq. ([l} . One may extract information about the gluon-decay stage 
in the MLLA region from gluon jet data in e"'"e~ collisions. In order to obtain the energy 
dependence of the function ]^ Gluon ^ directly experimental data for N^iuon (^Ej^t) 

in e+e~ annihilation [3]. One can see from Fig. [2] that ]\[Giuon jg constant at about 
Ejet < 1 GeV and it grows as a power of Ejet at higher energies. From the available 
e~^e~ collisions data shown in Fig. [21 we found that the energy-dependence of the mean 
charged particle multiplicity of gluon-jet can be approximately described by 

(7V^'"°") cx E^j^^, with (5 = 0.6 0.7 for Ej^t > 0.85 1 GeV. (3) 

It is essential to stress again that such behavior also follows from the theoretical estimates 
in the 3NL0 pQCD in the MLLA scheme [TT]. Now, using Eqs. (I2I3I) and assuming the 
typical energy of the gluon jet to be of the order of average saturation scale, we obtain. 



4 



A. H. Rezaeian Low-x Meeting 2011 



z 

V 







e e data 

* quark 

• gluon, OPAL (boost) 
O gluon, OPAL (qqg-qq) 

- - e" ' 




0.25 0.5 1 



E [GeV] 

g,q 



16 32 



Figure 2: The mean charged hadron muhiphcity of unbiased gluon Nj^ and quark 
jets in e^e~ annihilation, as a function of the jet energy. The plot is taken from [3]. 
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where s is the center-of-mass energy squared per nucleon pair and A is free parameter to 
be fixed with other experiments like DIS at HERA. In the above we assumed that the 
saturation scale for pp collisions is Qs < 1 GeV and for AA collisions we have > 1 
GeV. We take for the parameter S, the average value S — 0.65 from Eq. ([3]). In Eq. (jU, 
the average value of A w 0.22 in the the effective saturation scale for pp collisions can be 
obtained from kx factorization results given in Ref. |8 . Then, the power-law behavior 
given in Eq. ([5]) for A A collisions comes naturally without any extra freedom. The 
simple Eqs. (I4I5I) are enough to describe energy-dependence of all existing data on hadron 
multiplicity from RHIC to the LHC, see Fig. [3] (right). Indeed the full calculation [3] 
(shown by orange line in Fig. [3]) based on the factorization shown in Fig. [U agrees with 
the above argument based on the dimensionality. It is seen from Eqs. (I3I5P that the 
factorization accounts for the most energy-dependence of the hadron multiplicity s'^'^^, 
while the MLLA gluon-decay cascade (extracted from e+e^ annihilation data) brings 
only about s^ o^e extra effect. Notice that in the above we assumed that the atomic 
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Figure 3: Right: The energy behavior of charged particle pseudo-rapidity per participant 
pair for central AA and non-singlet difFractive pp collisions. The right panel plot is taken 
from Ref. [3]. Left: Pseudo-rapidity distribution of charged particles produced at RHIC 
and the LHC. The experimental data from the LHC are from [11[T]. 



number or A (size of nuclei) dependence of the saturation scale is factorizable from 
energy. This assumption can be simply justified from the observation that the centrality 
dependence of the multiplicity at the LHC is very similar to RHIC up to an overall scale, 
see Fig. [H This means that the energy and centrality or A dependence of multiplicity 
(and saturation scale) can be scaled out, in agreement with our assumption. Therefore, 
data seems to role out any extra energy dependence to be generated due to the geometry 
and the centrality of collisions. In Fig.[3l we compare our published prediction [3] for the 
rapidity distribution of charged hadrons coming from the full calculation with the recent 
CMS data [2] . It is seen in Fig. [3] that our predictions is in excellent agreement with 
data. Our prediction for dNAA/drj at midrapidity for — 5% Pb+Pb collisions at 5.5 
TeV is 1897 ± 133 [3]. We also show in Fig. [3] our prediction for the rapidity distribution 
of charged hadrons multiplicity for Pb-f-Pb collisions at 5.5 TeV j3]. 

In Fig. [4] (right), we show {2/Npar){dNAA/dri) at midrapidity (where Npar is the 
number of participant for a given centrality) has the scaling property at different energies. 
Notice that one can already observe similar scaling property at RHIC, namely dNAA/drj 
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Figure 4: Right: The scaled pseudo-rapidity density as a function of number of partic- 
ipant at midrapidity for A A colhsions at various energies. Left: The pseudo-rapidity 
distribution at RHIC 0.2 TeV at different centrahties. The plots are taken from Ref. [3]. 



at fixed energy but different centrahties falls into a single curve upto a normalization 
factor, see Fig. |4] (left). Both scaling properties shown in Fig. |4] can be easily understood 
within the CGC picture and follows from simple Eq. ([2]). We expect that the centrality- 
scaling for dNAA/drj at a fixed energy will be also valid at the LHC. We predict that 
{2/Npar){dNAA/ drf) for 5.5 TeV AA collisions at midrapidity to be about q gyj-g ge times 
bigger than the corresponding one at 2.76 TeV. 

In Fig. [5] (right) we show the experimental data from the STAR collaboration [T^] 
for {pt) I \J [dN'^ / drf) / St as a function of centrality where (pt) is the average transverse 
momentum and St is the overlap area between the colliding nuclei in the transverse 
plane. In our approach, we have (pt) / {dN/ drf) /St ~ -^^^ where n ^ ]\fGiuon 
Qs > 0.85 1 GeV corresponding to the excess of charged hadron production in the 
presence of jet-decay effects. It is seen from Fig. [5] that the ratio (pt)/ \/ (dN/ drj) / St 
at RHIC Au-|-Au collisions decreases for more central collisions and higher energies in 
accordance with our model and in contrast to the KLN type approach [13l [5] . We expect 
that the ratio (pt)/ \/ {dN/drj) / St will be further suppressed at the LHC compared to 
RHIC, see Fig. [5] Moreover, in the KLN type approaches [13] we have (pr) ^ x where 
X = {dN / drf) / St while in our approach we have (pr) ^ a;:0-264 ^^j. ^j^g (,g^gg t]\a± the 
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Figure 5: Right: The ratio {vt)I \J {dN^ / drf) / St at various centralities at RHIC. The 
data was constructed from three experimental measurements, the average transverse 
momentum {px) of 7r~, dN'^/dj] and St [H]- Left: Average transverse momenta as a 
function of [dN'^ / drf) / as at RHIC. The experimental data are from [T^. 



saturation scale is Qs > 0.85 4- 1 GeV. In Fig. [S] (left) we show the average transverse 
momenta as a function of ^ {dN / drj) / St ■ The STAR collaboration [12] has found that 
the experimental data for the charged pion at different centralities can be described by 
(pt) ~ Po + 0.072: where the constant po = 0.29 GeV was obtained from a fit and may 
be interpreted as primordial transverse momentum. In Fig. (SJ we also show that a fit 
driven by our approach prefers a smaller primordial transverse momentum of about pion 
mass pq ^ 0.14 GeV (or even pq ^ 0) and it reasonably describes the same data. Notice 
that at small multiplicity for very peripheral collisions at RHIC energy entire saturation 
formulation is questionable. 

To conclude: we extracted the energy-dependence of the gluon-jet decay cascade from 
e+e^ annihilation data. We showed that the energy- dependence of about s^ o^e ^]-^g 
gluon-decay cascade (when the average transverse momentum of the jet becomes about or 
bigger than 1 GeV) is exactly what explains the different power-law energy-dependence 
of hadron multiplicity in A A compared to pp collisions at the LHC. This effect is more 
important for AA collisions where the saturation scale is larger and gives rise to an extra 
contribution about 20 — 25% to the multiplicity in AA collisions at the LHC. On the 
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theory side, this emphasizes on the importance of higher order corrections beyond the 
leading log approximation in the kx factorization, and on general, this also gives rise to 
the outstanding problem that how the fragmentation processes can be accommodated 
within the CGC/satm'ation framework. 
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